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the  scatter  problem .to  the  filler  task  of  solving  a linear  system.  For 
cortaln  basic  topics  of  surface  profllos,  Including  the  sinusoidal 
profile  considered  here,  the  coefficients  of  the  linear  system  are: 
obtained  as  closed  fom  expressions  In  well-known  functions  (Bessel 
functions  fbr  slnusodlal  profiles  and  exponential  functions  for  piece- 
wise  linear  profiles).  The  theory  Is  thus  amenable  to  efficient  computer 
evaluation.^ 

In  part  I of  this  report  the  theory  Is  presented  In  detail;  part  II 
reports  on  numerical,  results  obtained  by  coiuter  evaluation  of  the  theory 
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SCATTERING  BY  PERIODIC  SURFACES  WITH  SINUSOIDAL 
HEIGHT  PROFILE-A  THEORETICAL  APPROACH 

PART  I : THEORY 
1.  INTRODUCTION 

Many  of  the  procedure*  used  to  study  scattering  from  periodic  surfaces 
are  baaed  on  the  original  approach  used  by  Lord  Rayleigh  in  1907  [1]  in  which 
a discrete  spectrum  of  outgoing  plane  waves  is  assumed  for  the  scattered 
field.  Difficulties  occur  in  applying  boundary  conditions  since  this  form 
of  the  scattered  field  does  not  necessarily  apply  at  the  scatter  surfaces 
[2,3J.  Other  approaches  to  study  scattering  from  periodic  surfaces  include 
perturbation  techniques  [4,5]  variational  methods  [6]  and  the  physical  optics 
or  Kirchhoff  approximation  [7,8,9].  All  of  these  methods  are  of  limited 
validity.  Only  recently, with  the  availability  of  high  speed  computers, has 
the  evaluation  of  rigorous  formulations  of  these  problems  become  practical  and 
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reliable  reeulta beenobtalned  [10-15].  A*,  excellent  survey  of  work  done  on 

•c***#r*Dt  ^roB  periodic  surfaces  can  be  found  in  the  report  by  Tong  and 
Senior  [13]. 


The  present  study  uses  an  independent  approach  to  the  problem  of  scatter- 
ing from  periodic  netal  surfaces.  It  is  basically  an  adaptation  and  exten- 
sion of  the  physical  optics  approxlaatlon  into  a rigorous  theory  which 
applies  to  surfaces  with  radii  of  curvature  not  necessarily  large  conpared 
to  wavelength.  The  specific  problea  being  considered  is  that  of  scattering 
o plane  waves  froa  a sinusoidally  varying  perfectly  conducting  surface.  The 
currant  Induced  in  this  surface  (which  produces  the  scatter  radiation)  la 
assuaed  to  have  the  fora 


K - K F 
P 


Kp  the  physical  optics  approximation  of  this  current  - or  a suit- 
ably chosen  aodification  of  it  - and  F is  a Fourier  Series  expansion  with 
a fundamental  period  equal  to  the  period  of  the  surface  (to  satisfy  Floquet's 
theorem  for  periodic  structures).  By  inposing  the  so  called  extended  boundary 
condition,  i.s. , the  physical  constraint  that  the  field  radiated  by  the  in- 
duced surface  currents  into  the  lower  (shielded)  half-space  cancels  the  inci- 
dent plane  wave  in  this  region,  a linear  system  is  obtained  for  the  complex 
Fourier  coefficients  of  the  current  distribution.  The  matrix  elements  of  this 
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system  and  the  Inhomogeneous  terms  are  Integrals  which  can  be  evaluated  lead- 
ing to  closed  form  expressions  in  terms  of  Bessel  functions.  Numerical  eval- 
uation of  the  theory  thus  Involves  solution  of  a linear  system.  But  since 
the  coefficients  of  this  system  are  closed  form  expressions rno  time  consum- 
ing integrations  are  involved  in  setting-up  the  coefficient  matrix.  Solu- 
tion of  this  system  (after  appropriate  truncation)  then  allows  the  determina- 
tion of  the  scattered  field  in  the  upper  (illuminated)  half— space  away  from 
the  conducting  surface;  the  complex  amplitudes  of  the  plane  wave  spectrum 

this  field  are  series  of  Bessel  functions.  By  relaxing  the  degree  of 
truncation  of  the  linear  system  in  the  numerical  evaluation  (i.e.,  by 
Increasing  the  order  of  the  system  actually  evaluated) ,it  should  in  principle 
be  possible  to  arbitrarily  Increase  the  accuracy  of  the  approach,  although  this 
has  not  yet  been  demonstrated. 

One  of  the  reasons  for  conducting  the  present  study  was  to  Investigate 
whether  or  not  the  specular  reflection  coefficient  of  periodically  corrugated 
metal  surfaces  is  significantly  dependent  on  the  direction  of  polarization  of 
the  incident  field.  Thus  both  TE—  and  TM— polarization  of  this  field  are 
considered  where /TE-polarization  denotes  the  case  that  the  electric  field 
strength  is  directed  parallel  to  the  direction  of  the  surface  grooves  and 
TM-polarization  the  case  that  the  magnetic  field  strength  has  this  direction. 
Interest  in  this  polarization  question  stemmed  from  the  need  to  verify  the 
superiority  of  TM-  over  TE-polarization  in  minimizing  false  guidance  of  micro- 
wave  scanning  beam  landing  systems  [16],  Reflections  from  large  metal 
structures  near  runways  such  as  hangars  (which  usually  have  periodically 
corrugated  walls  and  doors)  or  rows  of  parked  airplanes  are  a likely  cause  of 
misinformation  received  by  landing  aircraft,  and  the  question  arises  if  such 
reflections  can  be  reduced  by  appropriate  choice  of  polarization.  The 
numerical  results  confirm  experimental  evidence  that  TM-polarization  in 
general  leads  to  substantially  less  specular  reflection  than  TE-polarization,  in 
particular  in  the  practically  interesting  range  of  low  incidence  angles  (near 
grazing) . At  grazing  incidence  itself  the  specular  amplitude  reflection  co- 
efficient of  course  is  -1,  for  both  TE-  and  TM-polarization. 

Por  a summary  of  the  new  theory  we  refer  to  a preceding  shorter  report 
[17]  which  also  discusses  the  case  of  circular  polarization  of  the  incident 
field.  It  is  shown  In  this  report  that  circular  polarization  is  highly  effec- 
tive in  reducing  higher  order  grazing  lobes  with  the  same  polarization  (cir- 


[16]  Demko,  P.  S.,  (1972),  "Polarization/Multipath  Study.",  Tech.  Memorandum 
VL-5-72 , Avionics  Laboratory , U.S.  Army  Electronics  Command  Ft.  Monmouth, N.J. 

[17]  Schwering,  P.  and  Whitman,  C.,(1977)  "Scattering  by  Sinusoidal  Surfaces", 
Tech.  Kept.  E COM- 4496,  Communications/ADP  Laboratory,  U.S.  Army  Electronics 
Command,  Port  Monmouth,  N.J. 
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cular  with  th«  same  sense  of  rotation)  es  the  Incident  wave.  On  the  other 
hand,  grating  lobes  with  the  opposite  sense  of  rotation  are,  in  general, 
present  with  significant  amplitudes.  Such  grating  lobes,  however,  would  not 
be  received  by  an  airborne  antenna  polarised  for  optimum  reception  of  the 
primary  microwave  beam. 

2.  FORMULATION  OF  THE  PROBLEM 


Consider  a perfectly  conducting  sinusoidally  varying  surface  defined  by 
the  relation 

x 

z - f(x  ) - h sin  (2w  — r-  ),  - • < x , y < •. 

O O u O O 


A uniform  plane  wave  Is  assumed  Incident  upon  this  surface  at  an  angle  of 
9 degrees  with  a suppressed  time  dependence  exp  (lent) ; see  Fig.  1.  The 
field  strength  components  of  the  primary  wave  are  In  the  case  of  TE-polar- 
isatlon 

_p  ik(x  sine  - s cos0)  (2a) 

* 6 


and  In  the  case  of  TM-polarlsatlon 


_p  -ik(x  sin 9-  z cos0) 
H ■ 6 


ik  3x 


(2b) 


In  the  TR-case,the  electric  field  strength  of  the  primary  wave  Is  assumed  to 
have  unit  amplitude  and  in  the  TM-case, the  magnetic  field  strength;  k a Il- 
ls the  wave  number.  * 


Aa  it  Is  well  known  from  the  theory  of  optical  gratings,  the  Incidence  power 
Is  scattered  into  a finite  number  of  discrete  directions  6^ 

sinO  ■ slnfl  + m -4  (3a) 

n a 

where  X is  the  wavelength  and  m is  any  integer  In  the  range 

- ^(1  + sine)  < m < y(l  - sln0)  (3b) 

The  scatter  f leld^ In  other  words, comprises  a finite  spectrum  of  propagating 
plane  waves  carrying  the  real  power  reflected  by  the  surface.  This  spectrum 
Is  supplemented  by  an  Infinite  but  discrete  set  of  evanescent  plane  waves  as- 
sociated with  the  Integer  m values  outside  the  range  (3b).  For  such  m- 
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value*  the  propagation  angles  (3a)  become  complex,  describ 
waves  exponentially  decreasing  In  the  positive  z-direction 
TE-polarlxation, the  scatter  field  thus  takes  the  form 


and  in  the  case  of  TH-polarization 


(1)  ~ik(x  sin0  + z cos0  ) 


"n  5 6 “"Plitudes  of  the  component  plane  waves 

Using  (3a)  it  can  be  seen  that  these  component  waves, 
actor  . itaains,  periodic  In  x with  the  period  d of 
Floquet  s theorem).  They  have  therefore  been  termed 


3.  INTEGRAL  REPRESENTATION  OF  SCATTER  FIELD  AMPLITUDES 

In  the  case  of  TE-polarization,the  incident  field  h 
currents  in  the  metal  surface.  The  electric  field  strei 
field  can  be  derived  from  this  (unknown)  current  distribi 
to  the  relation 


(kr)dS, 


where  3/3n  denote*  the  derivative  normal  to  the  surface  S_. 

O u 

Since  the  Integrals  in  eqs.  (Sa)  and  (5b)  extend  over  an  infinite  range, 
they  ere  not  emanebl*  to  numerical  evaluation.  However,  because  of  the 
periodicity  of  the  surface  profile,  these  integrals  can  be  transformed  into 
integrals  extending  over  one  period  only.  Using  Floquet's  theorem  we  conclude 
that  the  surface  current  densities  KyOO  “»d  Kt(xQ)  are  periodic  with  d, 
apart  from  a phase  factor  exp (-ikxosin07 determined  b-r  the  incident  plane 
wave . Hence 


where 


VV  ‘ 


-ikx  sin8l  , . 

e ° Ky( 


v , . -iVoc  Sinflj:  , * 

Kt(xo)  - e o Kt(xo) 


KyCx^)  - KyW) 
Kt(xo+d)  - Kt(d) 


(ba) 

(6b) 

(6c) 

(6d) 


Subdividing  Integrals  (5a)  and  (5b)  into  sums  of  integrals,  each  extending 
over  one  period  only  of  the  scatter  profile,  we  obtain  after  interchanging 
summation  and  integration 


E*(x,z) 


-d/2 

- r J 

-d/2 


z-z)K(x) 
o y o 


H*(x,z) 


i 

4 


z-z  ) K (x  )dx 
o c o o 


(7a) 

(7b) 


where  g^  and  g^  are  given  by 


gB(x-xo,*-zo)  - £ Ho<2)  [k^(»-x0-nd)2  + (*-*0)2)  ]}e" 
n— »> 

• — »>  -{(■&  hr  k)Z_  H»2>tk  { 

+ («-.0>2  j" 


lnkd  sinO 


(7c) 


(7d) 


-lnkd  sinO 
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We  have  expressed  the  surface  element  ds  and  the  differential  3/3n  (in 

the  expression  for  g„  in  terms  of  x and°z  : ° 

n o o 


ds. 


■ 


dx 


9n 


vr. 


— (z*  -JL- § \ 

,,2Uo  3x  3e> 

z o o 

o 


where  z and  z are  the  height  and  slope  respectively  of  the  surface  profile. 
Hence,  with  eq°(l) . 


z ' ■ t ° ■ 2*  cos  (2*  ~ ) 

o dx.  d d 


(8) 


We  now  apply  Poisson's  summation  formula  [18],  [19]  which  we  write  in  a 
form  appropriate  for  our  present  purpose: 


J H^2>[k  {(x-xo-md)2  + (z-zo>2  }ift-imkdalne 


(9) 


- Ycoi*  e‘lk((Mo)sln^,-,olcoseJ 
kd  ^ * 


with 

sin0  ■ sin6  + m X/d 

n 

(9a) 

and 

Re  (cos6^>0,  Im(cos0a)<O 

(9b) 

Then  equations  (7a)  to  (7d)  reduce-after  Interchanging  integration  and  sum- 
mation- a second  time-te  space  harmonics  representations  in  accordance  with 
•V*  (^*)  (♦*>).  In  the  case  of  TE-polarization  we  obtain 


. V*  »(1)  -ik(x  si°Q  + * cos0  ) (10a) 

ay  4_,  a e * m 


[18]  Pel sen,  L.  B.  and  Marcuvits,  N.  (1973),  "Radiation  and  Scattering  of 
Waves",  Prentice-Hall,  Hew  Jersey. 

[19]  Morse,  P.  M. , and  Feshbach,  H.  (1953),  "Methods  of  Theoretical  Physics- 
Part  I,"  McGraw-Hill,  New  York. 
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i f ° A r if  r v ik(x  sinO  + z cos6  ) f 5 

2d  V e ^80  y(xo)e  0 ■ 0 »Vl+z’2  dx 

o m j 7 o ( 

J-d/2 

•nd  in  the  case  of  TM-polarization 

00 

H8  . V o-ik(X  8in®  + Z CO80  ) 

y Z_»  o m m 


Wlth  n d/2 

a ^ f if  /„  \ ik(x  sin Q + z cos0  ) f,  * 
m 2d  I Kt  o)e  0 « o ® i1_zo  t®n0m]dxo 

-'-d/2 

In  both  cases  we  assume  z * h - Max  (zq)  to  insure  convergence,  and 

Re(cos0  ) > 0,  Im(cos0  ) < 0 
in  m 

to  satisfy  the  radiation  condition  (see  eq.  (?b)). 


o 

(11a) 

(10b) 


(lib) 


Note  that  while  eqs.  (7a)  to  (7d)  represent  the  scatter  field  at  all 
points  above  the  scatter  surface,  including  the  space  range  within  the 
grooves,  the  space  harmonics  representation  (10a)  and  (10b)  converge  in  the 
space  range  z z h but  may  diverge  within  the  grooves.  Hence  this  repre- 
sentation cannot,  in  general,  be  used  to  formulate  the  boundary  condition 
Etang  " 0 at  the  metal  surface  (Rayleigh  Hypothesis  [20A  For  our  present 
purpose  the  space  harmonics  method  is  appropriate:  the  scatter  problem  is 
solved  using  an  approach  which  does  not  require  explicite  formulation  of 
this  conventional  boundary  condition. 


4.  CALCULATION  OF  INDUCED  SURFACE  CURRENT  DENSITIES 

Equations  (10a)  to  (lib)  provide  a representation  of  the  scatter  field 
of  the  periodic  surface  in  terms  of  the  (unknown)  current  distribution  induced 
incldent  plane  wave‘  To  determine  this  current  distribution  we  use  the 
following  approach.  The  metal  surface  is  assumed  removed  while  the  current 
distribution  existing  in  this  surface  is  maintained  in  place.  If  an  ap- 
proximation is  used  for  this  current  distribution,  deviating  from  the  actual 
current  density,  then  a field  will  exist  not  only  in  the  half-space  above  but 
also  in  the  half-space  below  the  current  carrying  surface.  This  latter  field 
can  be  derived  by  a procedure  similar  to  that  employed  in  calculating  the 
field  in  the  upper  half-space;  the  condition  that  this  field  must  be  identi- 
cally r®ro  then  utilized  to  determine  the  correct  current  distribution 
on  the  metallic  surface. 


Equations  (7a)  to  (7d)  hold  for  z > 0 as  well  as  for  a < 0.  Inserting 


[20]  For  a discussion  of  Rayleigh's  method  and  it's  limitations  see  Tong  and 


Poissons  summation  formula  (9),  while  observing  that  In  the  lower  half-space 
* K *0*  we  obtain  In  the  case  of  TE-polarization 


PT  p(2)  -ik(x  sine  - z cos8  ) 

y Zj  m e m m 


t<2>  - - 


_i_  nr _i_  rdn, 

2d  V e0  coaea  I lty<’ 
-d/2 


ik(x 

e o 


sin9  - z cos6  ) 
mo  m 


(12a) 


1+z ,2  dx 
o o 


(13a) 


and  in  the  case  of  TM-polarization 


«v  -y 

y A m 


(2)  -ik(x  sin8  - z cos6  ) 


(12b) 


H«>  - - 


-JJ  f ■tt(«0)elk<Xo  ’ln9»  ' 'o  “‘VlW  t.n8m)d» 

^ ~d/2  Q3b 


(13b) 

As  before,  the  propagation  angles  6m  are  given  by  eqs.  (9a)  and  (9b).  To 
Insure  convergence  It  is  assumed  that  z $-h  » Min  (z  );  for  z and  z’  see 
eqs.  (1)  and  (8) , respectively.  In  the  lower  as  in  ?he  upper°half-s$ace  the 
scatter  field  is  thus  obtained  as  a superposition  of  space  harmonics  i.e., 
as  a discrete  spatial  spectrum  of  propagating  and  evanescent  plane  waves 
whose  x-dependence-except  for  a common  phase  facter-ls  periodic  with  the 
period  d of  the  surface  profile.  The  propagation  directions  0 of  the  upper 
and  lower  half-spaces  correspond  to  each  other;  they  are  related  by  imaging 
at  the  plane  z • 0. 

The  current  distribution  existing  on  a metal  surface  is  now  obtained  from 
the  condition  that  the  field  in  the  lower  half-space  must  be  identically  zero. 
This  means  that  the  space  harmonic  ofi  order  zero  must  cancel  the  incident 
plane  wave  (it  travels  in  the  same  direction)  and  that  all  higher  order  space 
harmonics  must  vanish.  Hence  for  TE-polarizatlon 


5<2>  - 


and  for  TM  - polarization 


0 for  m ^ 0 

1 for  m ■ 0 


(14a) 


0 for  m + 0 

-1  for  m ■ 0 


(14b) 


Th«a«  aquations  in  conjunction  with  eoa  /nv>  * 

«us=i  szxr.  ajLSL. 
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cos0 


fj  a"ik(x0  alne  ' *0  CO80) 
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Ct(xQ)  - - 2 e~ik(xo  flln0  ' *o  cos0^Ft(xo) 


(15a) 


(15b) 


* f'md**“t*1  1-1  » th. 
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(15c) 
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R _ n " f *y  ol>t*lnad  by  formulating  tha  boundary  condition 

~ * rS&s'.ra:  aragf 


121] 


Sr*?CiSrdujfofe®thSPS25i2fi exP^®88ion  for  Ky 

ai settle  fiSfd  •tr•1^8t805?af^2a5^?Sa^^iveU?Attf?§P6t, 
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■•gnatic  fiald  strength  in  tha  TM-case. 


but  not  in  chat 
of  the 

cnee  and  of  the 


[22]  Th«  physical  optics  approximation  would  raault  if  the  factor 

function  (phaae  tarn)  in  <**>  — «- 


[23]  See  reference  [13]  0n  page  2. 
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Inserting  eqs.  (15a) , (15b)  and  (15c)  into  eqs. -(13a)  and  (13b)  and 
performing  the  integrations  leads  to  closed  form  expressions  in  terms  of 
Bessel  functions.  Thus, 


CE  ei2ir[(m+n)  + — sin(2ir-j)  (cosB-cosO^)] 


cosfl 

cosG 


n-fm+1 


8in(2ir-j)  (cos0-cos0  )] 


2V  m+n-t-l 


with 


where  the  angles  0 are  defined  by  eqs.  (9a)  and  (9b).  With  these  expressions 
for  and  H'^),  conditions  (14a)  and  (14b)  reduce  to  the  following 

linear  systemsmfor  the  unknown  current  coefficients  CT  and  C*J: 


Bqa.  (17b)  and  (17d)  yield  explicit  expressions  for  the  aero  order  coeffi- 
cients C?  and  C , reapectively : 


Ualng  these  expressions  in  eqs.  (17a)  and  (17c),  the  linear  systems  for  the 
higher  order  current  coefficients  take  the  formt 


Solution  of  the  linear  systems  (18c)  and  (18d)  requires  a matrix  inver- 
sion which  is  efficiently  handled  by  available  computer  routines.  A parti- 
cular advantage  is  that  all  matrix  coefficients  and  Inhomogeneous  terms  are 
closed  form  expressions  in  well-known  functions  so  that  no  time  consuming 
numerical  integrations  are  involved  in  computing  the  coefficients  of  the 
matrix  to  be  Inverted,  though  a routine  for  Bessel  functions  of  complex 
arguments  is  needed.  Inversion  of  course  requires  truncation  of  the  in 
principle  infinite  linear  systems.  We  assume  without  proof  that  the  current 
distributions  computed  will  converge  towards  the  actual  distributions  when 
the  degree  of  truncation  la  relaxed,  i.e.,  when  the  order  of  the  linear 
systems  is  step  by  step  Increased.  As  the  critique  of  the  Rayleigh  method 
has  shown  such  assumptions  have  to  be  regarded  with  caution.  In  contrast  to 
the  Rayleigh  hypothesis,  there  are  no  obvious  physical  reasons  for  doubting 
the  validity  of  the  present  approach;  but  problems  have  been  encountered  in  the 
mmwrlcel  evaluation  when  the  groove  depth  (2h)  exceeded  1 to  2 wavelengths. 


(See  part  II  of  this  report) . 


5.  SERIES  REPRESENTATION  OF  SCATTER  FIELD  AMPLITUDES 


EH 

When  the  Fourier  coefficients  C_  and  C , of  the  current  distributions 


v_  WMV  v I Wfc  WUiiCUb  UXOLiXUUUlUUD 

Ky  and  Kt  are  known,  evaluation  of  tne  plane  wave  amplitudes  of  the  scatter 
field  in  the  upper  half-space  is  straight  forward.  In  the  case  of  TE-polar- 
ization  we  obtain  by  insertings  eqs.  (15a)  and  (15c)  into  eq.  (11a) 


E(l)  „ _ 1 cos 6 fd/2  y1  CE  ei2Tr|m+ii)^2.  + ^ sln(2w  ~)  (cos0  + cos6ml]  dx 
m d cos0mJ_d/2  L->  « o 


^Trfcm+n)—  + y sln(2w  ~)  (cose  + cos6m1] 


(19a) 


>lm>C 


and  in  the  case  of  TM-polarization  by  using  eqs.  (15b)  and  (15c)  in  eq. 
(lib) 


d/  2 

[j(l)  i T ei2irC(®+n)  -jp  + y sin(27r  -p)(cos0  + cos0  )] 

m d J-d/2  n m 


(19b) 


*(l-2ir  | cos(2Y)  tan9m)  dx0 


!<■> 


.m+n+1 


[1  + W°i»>  c" 


(19c) 


where  in  both  cases 


t. 

alm  " y(c°80+  CO80m)  ^ Pm  = 2ir  tan0ffl 


(19c) 


We  have  used  here  eqs.  (1)  and  (8)  for  zQ  and  z'  , eq.  (9a)  for  sinQm,  and 
have  furthermore  assumed  that  the  Fourier  series  of  the  current  distribu- 
tions converge  sufficiently  to  permit  interchanging  summation  and  integra- 
tion. 


6.  CONSERVATION  OF  POWER  DENSITY  AND  RECIPROCITY  THEOREM 


Two  physical  criteria  are  used  to  check  the  validity  of  our  numerical 
results.  They  are  conservation  of  power  and  reciprocity.  The  former  requires 
that  the  power  reflected  per  unit  area  from  a perfectly  conducting  surface 
equal  the  power  incident  per  unit  area  upon  the  same  surface,  i.e., 


where 


-Inc  _refl 

sj-  »j. 

sx-  J Rs(E  x H*-n) 

is  the  normal  component  of  the  Poynting  vector  and  n li  a unit  vector  per- 
pendicular to  the  surface.  In  the  present  case  it  is  convenient  to  formu- 
late this  condition  for  one  period  of  the  metal  surface,  utilising  a plane 
horlionttl  test  surface  subtended  between  two  adjacent  maxima  of  the  surface 
profile  (to  insure  convergence  of  the  space  harmonics  representation  of  the 
scatter  field).  Por  the  primary  power  incident  through  the  test  surface  we 
obtain  with  eqs.  (2a)  and  (2b) 


for  TE-polarisation  (20a) 

for  TM-polarisation 

and  for  the  reflected  power  propagating  through  this  surface  we  obtain  by 
using  eqs.  (4a)  and  (4b) 


for  TE-polarisation  (20b) 


for  TM-polarisation 


where  the  summation  is  limited  to  the  propagating  space  harmonics  i.e.,  to 
spectral  orders  m satisfying  condition  (3b).  Eq.  (20b)  expresses  the  fact 
that  the  space  harmonics  are  mutually  orthogonal  over  one  period.  Conserva- 
tion of  power  thus  requires 


cos6  «•  < 


|K 

En 


(1)  1 2 


(1)  1 2 


cose. 


cose. 


for  TE-polarisation  (20c) 


for  TM-polarisation  (20d) 


providing  a criterion  for  checking  the  accuracy  of  the  computed  results. 


The  second  physical  condition  tested  is  reciprocity.  As  shown  in 
appendix  A,  the  Lorents  reciprocity  theorem  applied  to  the  problem  of 
periodic  surface  scattering  may  be  stated  as  follows:  If  a plane  wave  inci- 
dent from  direction  0®  produces  a space  harmonic  of  order  m propagating  in 
the  direction  ©£  • - 6 , then  a plane  wave  arriving  from  direction  0b  pro- 
duces a space  harmonic  of  the  same  order  m traveling  in  direction  - -0* 
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[24]  and  the  complex  amplitudes  of  the  two  space  harmonics  are  related  by 


CO80 


for  TE-po lari sat Ion  and  by 


(Hm>  -se;-  (H-'rcose°  (21b) 

for  TM-polarization.  We  have  assumed  here  that  the  two  Incident  waves  have 
equal  amplitude  and  phase. 

The  power  and  reciprocity  criteria  supplement  each  other.  While  the 
former  criterion  checks  on  the  accumulated  power  of  all  propagating  space 
harmonica,  the  reciprocity  criterion  tests  the  amplitude  and  phase  of 
individual  space  harmonics. 


APPENDIX  A 


Reciprocity  Relation  for  Periodic  Surface  Scattering 

Two  TE-polariied  plane  waves,  with  propagation  angles  0a  and  6 


ik(x  sin0  -z  cos0a) 


£ -ik(x  sinO  -z  cos©*) 


cos6  + e 


ik(x  8in0  -zcos0D) 


£ -ik(x  sln0  -zcos0°) 


cos0  + e 


124]  The  minus  signs  In  the  relations  ©£  - -e"  and  0^  - - 6a  a 
definition.  Propagation  angles  of  Incident  and  scattere 
going  and  outgoing  waves)  are  counted  with  opposite  signs 


assumed  incident  on  a periodic  aetal  surface  vith  height  profile  (1) 
i.e. , 


The  respective  scatter  fields  are  denoted  by  K , H and  B , H°  . In 
the  space  range  above  the  surface  grooves,  * > Sax  ^ 0 ) * h?Cthese  scatter 
fields  can  be  expreesed  by  their  apace  harmonic  representations: 


,a  -ik(x  sine*  + xcosg*) 


+ icoaO  ) 


where  e , e , e are  the  unit  vectors  of  the  x,y,«  -directions,  respec 
tlvely,  and”from  eq.  (9a) 


The  Lorenta  reciprocity  theorem  applied  to  the  tvo  dimensional  problem 
considered  have  states  that 


where  E , H and  E , H are  arbitrary  electromagnetic  fields,  S is  any 
closed  contour  not  containing  (active)  sourcee,  and  n la  a unit  vector 
normal  to  S.  We  chose 


and  apply  the  reciprocity  relation  to  the  contour 


shown  in  Fig.  A. 2 where  has  the  length  d of  one  period  of  the  metal  sur- 
face, and  S3  have  arbitrary  height,  and  S4  follows  closely  the  metal  sur- 
face. Clearly,  the  contribution  of  Sa  to  the  reciprocity  integral  RI  is 
rero  since  at  the  metal  surfr.ee  Ea,  E&  - 0.  On  the  remaining  portion  of 
the  contour  the  scatter  fields  may  be  expressed  by  their  space  harmonics 
representations  (A. 4)  and  (A. 5). 


We  now  make  the  assumption  that  the  direction  of  arrival  of  the  inci- 
dent field  Eb,  coincides  with  the  propagation  direction  of  one  of  the 

space  harmonics  (spectral  order  m - m')Df  the  scatter  field  Ea  , Ha  . In 
other  words,  {25]  ”8C  "®c 


-0b  - 0a, 

m 


(A.  8) 


Then,  as  can  be  shown  from  eqs.  (A. 6),  the  propagation  direction  of  one  of 

the  space  harmonics  of  the  scatters  field  E^  , H*5  will  by  necessity 

coincide  with  the  direction  of  arrival  of  iu8idenlcf ield  Ea  , H*  : 

-inc  -inc 


-9*  - 9b, 
m 


(A. 9) 


The  order  of  this  space  harmonic,  m - m',  is  the  same  as  in  eq.  (A. 8)  and 
from  (A. 6)  and  (A. 8) 


- Y (sin9a  + sin0b) 


(A. 10) 


With  a later  application  in  mind  we  note  that  if  eqs.  (A. 8)  and  (A. 9)  are 
satisfied,  the  propagation  directions  of  scatter  field  (b)  are  those  of 
scatter  field  (a)  imaged  about  the  z - axis,  and  vice  versa.  This  statement 
is  verified  by  adding  eqs.  (A. 6)  and  utilizing  eq.  (A. 10)  to  show  that 


sin0a  + sin0b  - (m+n-m' ) ~ 
m n d 


from  where  it  follows  that 


0b,  - - 0a 

m -m  m 


(A. 11) 


(A. 12) 


Using  eqs.  (A.l),  (A. 2),  (A. 4)  and  (A. 5)  we  now  formulate  the  reciprocity 
Integral  (A. 7): 


1251  With  regard  to  the  minus  signs  in  eqs.  (A. 8)  and  (A*9) 
angles  of  arrival  (0a,0b)  and  scatter  angles  (0j,  ej[) 
opposite  signs. 


see  Fig.  1; 
are  counted  with 
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/ • m n 


[-e(sin0a-sin0b)-e  (cos0a-cos6b) 
—x  in  n -x  m n 


•ndS 


(A. 13) 

Ha  have  utilised  here  eqs.  (A. 6),  (A. 10)  and  (A. 11)  to  write  the  exponents 
of  the  integrand  in  a fore  which  shows  that  the  Integrand  is  periodic  in  x 
with  the  period  d of  the  Betel  surface.  As  a consequence  of  this  periodi- 
city, the  contributions  of  and  to  RI  cancel  each  other  (observe 
that  n «-ex  on  S.  and  n • ♦ e on  S^).  Furthermore,  only  those  terns  of 
the  integrand  whose  x-dependent  phase  factors  reduce  to  unity  yield  a non- 
sero  contribution  when  integrated  along  S2. 

Thus: 


RI 


• (cos0fl 

The  last  tern  within  the  parenthesis  (sun  tern)  vanishes  because  of  eq.  (A. 12). 
The  reciprocity  theoren  (A.7)  thus  reduces  to  the  condition 


E*,  cos0a,  • Eb,  cos0b, 


(A. 15) 


This  result  nay  be  stated  in  the  following  fashion.  If  a plane  wave  Incident 
froa  direction  0*  produces  a space  harmonic  traveling  in  direction  -6**,  then 
a plane  wave  incident  fora  direction  0“  will  generate  a space  harmonic  of  the 
seme  order  propagating  in  direction  -6a.  If  the  two  incident  waves  moreover 
have  the  same  amplitude  and  phase  then  the  two  space  harmonics  considered 
also  have  equal  phases  while  the  ratio  of  their  amplitudes  is  equal  to  the 


inverse  ratio  of  their  direction  cosines 


Application  of  the  reciprocity  theorem  to  TM-polarized  fields  follows 
same  pattern  as  used  here  for  the  case  of  TE-polarization.  The  result 


corresponds  in  all  respects  to  condition  (A. 15).  H®,  and  H , are  the 
complex  amplitudes  of  space  harmonics  propagating  in  directions  -0m  = 
and  -0D  = 03,  respectively,  generated  by  incident  plane  waves  of  unit 
plitude,  zero  phase  and  respective  angles  of  arrival  0a  and  0. 

APPENDIX  B 


Summary  of  Formulas 


With  regard  to  the  coordinates  and  surface  parameters  we  refer  to  Fig. 
of  the  main  text.  If  the  groove  width  of  the  periodic  surface  is  suffi- 
ciently large,  d » A,  then  it  is  reasonable  to  assume  that  the  current 
density  in  this  surface  is  primarily  determined  by  local  effects.  In  other 
words,  it  may  be  assumed  that  the  current  in  each  surface  element  is  the 
same  as  if  this  surface  element  were  part  of  an  infinite  metal  plane 
(tangent  plane)  having  the  same  slope  as  the  surface  element.  Then  for 
TE-polarization 


and  for  TM-polarization 


ik(xD  sin0-zo  cos0) 


where  1^,  Hy,  evaluated  at  x0,z0  are the  magnetic  field  strength  componei 
of  the  incident  plane  wave  at  the  metal  surface,  0 is  the  angle  of  incidei 
and  zD  and  z'Q  are  the  height  and  slope  of  the  surface  profile  as  given  by 
eqs.  (1)  and  (8): 


h sin(2ir 


In  the  TE-case  the  electric  field  strength  and  in  the  TM-case  the  magnetic 
field  strength  of  the  incident  plane  waves  have  unit  amplitudes. 


Th*  co«pl«x  amplitudes  of  the  apace  harmonica  of  the  scatter  field  are 
obtained  by  using  current  distributions  (B.l)  and  (JJ.2)  In  eqs.  (Ha)  and 
(lib)  respectively.  Thus  for  TE-polarlzatlon 


2 xo 

KC1).  _ 1 coab  [ l(2wm  --  + kz0(cosd  + cose-)...  , . . 

m d cose.  J e d ■ (l+*o  tan0)  dx° 


(B.  3) 


and  for  TM-polarlzatlon 


-I  / 3(2irm  + k*o  (co*9  + cose.),  , . x . 

d 1 d " (l-z£  tan 0m)  dxo 

n 


(B.  4) 


where  6a  la  defined  by  eq.  (9a).  If  shadow  effects  are  taken  Into  account, 
the  lnte*ratlon  la  limited  to  the  aubrange  In  which  the  metal  surface  is 

ahI^1Lii1UM^Uft*d*by  tH*  lncldent  Plane  wave  while  in  the  geometrical 
shadow  sonea  the  surfece  currents  are  assumed  to  be  zero:  See  Fie  B 1 

•*  l»«.r«Ion  1„  Chi.  c...  « 


XI  + h tan9  sin (2*  ^J)  - h\/tan2e  - (^-  £)*  + (x2  - d) 


*2  • I j1-  i co*'1(ir  i "89>| 


(B.  5) 


If  shadows  effects  are  neglected  or  no  geometrical  shadowing  occurs 
(l.e.»  for  0 < tan’H^  y) ) the  Integration  la  performed  over  a full  period 
of  the  surface  profile  i.e.,  over  the  range  -d/2  * x*,  (d/2.  In  this  case 
the  Integrals  (B.3)  and  (B.4)  can  be  evaluated  In  closed  form  leading  to 
Identical  expressions  for  TE-and  TM-po lari ration: 


' Hl1>-  <-»tl  + <«1»> 


(».6> 


with  P, 


- (-1)* 


2ir  -r  tan0B 


1 + cos  (6  4-  9, 


(cos6  + cos0n)co86ra 


aia  ■ 2ir  ^(cosO  + cosen) 
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